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F.-J. Jiang 1 '!] and B. C. Tiburzi 2 '0 

1 Institute for Theoretical Physics, Bern University, Sidlerstrasse 5, CH-3012 Bern, Switzerland 
Maryland Center for Fundamental Physics, Department of Physics, 
University of Maryland, College Park, MD 20742-4111, USA 

We use a generalization of chiral perturbation theory to account for the effects of flavor twisted 
boundary conditions in the Breit frame. The relevant framework for two light flavors is an SU(6\4) 
partially quenched theory, where the extra valence quarks differ only by their boundary conditions. 
Focusing on the pion electromagnetic form factor, finite volume corrections are calculated at next- 
to-leading order in the chiral expansion and are estimated to be small on current lattices. 

PACS numbers: 12.38.Gc,12.39.Fe 



Twisted Boundary Conditions and the Breit 
Frame. Simulations of QCD on Euclidean spacetime 
lattices are making progress towards quantitative un- 
derstanding of strong interactions A source of sys- 
tematic error in these calculations is the finite lattice 
volume. While observables generally depend upon the 
lattice volume, there is a potentially more serious ef- 
fect: periodic boundary conditions on quark fields limit 
the available momentum modes to integer multiples of 
2tt/L, where L is the size of the lattice. This presents 
difficulty for the study of hadronic properties at small 
momentum. Periodic boundary conditions, however, 
are often chosen for convenience; and, a more general 
class of boundary conditions, twisted boundary condi- 
tions (TwBCs), see e.g. [2j, are possible. A TwBC by 
an arbitrary twist angle 9 on the matter field tp, has the 
form: tp{x + L) = e l9 ip(x), and consequently the matter 
field has kinematic momentum p — (2irn + 9)/ L, where n 
is an integer. The ability to produce continuous hadronic 
momentum has motivated the studyof TwBCs in lattice 
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Producing continuous hadronic momenta via TwBCs 
does not necessarily give one the ability to probe ampli- 
tudes at continuous momentum transfer, fn fact, these 
techniques are of no avail for flavor singlet form factors 
which require operator self-contractions. In a current 
self-contraction, the boundary condition chosen for the 
quark in the current does not affect the momenta of the 
external states. Flavor non-singlet form factors, how- 
ever, can be accessed at continuous momentum transfer. 
A conceptually clear way to see this is to consider form 
factors of flavor changing currents Q. If the current in- 
duces a change from a quark flavor satisfying periodic 
boundary conditions to one satisfying TwBCs, continu- 
ous momentum transfer, 6/L, will be induced. There 
is another way to utilize TwBCs for flavor non-singlet 
currents [Tij ]. For this method, imagine that the current 
strikes a quark of a given twist angle 9, and then produces 
a different twist angle 9'. Here the momentum transfer 
(9' — 9)/L is induced. In the infinite volume limit, these 
boundary conditions become irrelevant and the current 
thus produces momentum transfer by striking a single 



(non self-contracted) quark. Furthermore one can choose 
9' = —9 to work in the Breit frame, which often simpli- 
fies the calculation of form factors. At finite volume, 
however, the initial- and final-state quarks are actually 
distinct because a given field can only have one boundary 
condition. Hence the lattice action must be described by 
a theory with an enlarged valence flavor group [l4] . The 
purpose of this note is to utilize the effective field theory 
for the Breit frame implementation of TwBCs, discussed 
in [l4j |. to compute volume corrections to observables, in 
particular the pion electromagnetic form factor. 
Partially Quenched Chiral Perturbation Theory. 
To address the long-distance effects of TwBCs in the 
Breit frame, we must formulate the low-energy effective 
theory. From the discussion in 14], this variant of chi- 
ral perturbation theory includes additional fictitious va- 
lence quarks. These fictitious quarks only differ by their 
boundary conditions. We can remove the twists by a field 
redefinition in favor of periodic quark fields. In terms of 
these fields, the quark part of the partially twisted, par- 
tially quenched QCD Lagrangian appears as 



C = Q[lp+m Q ) Q 



(1) 



where the vector Q accommodates ten quark fields, 
Q = (ui,U2,di,d2,j,l,ui,U2,di,d 2 ) T . The field Q 
transforms under the fundamental representation of the 
graded group SU(6\4). The quark mass matrix in the 
isospin limit 1 of the valence and sea sectors is given by 
niQ = dia,g(m u ,m u ,m u ,m u ,mj,mj,m u ,m u ,m u ,m u ). 
Because we implement twisted boundary conditions only 
in the valence sector, we have additionally chosen to work 
away from the unitary point, where rrij = m u . With this 
choice, we can track sea-quark effects with ease and ac- 
count for partial quenching errors. Formulae relevant at 
the unitary point can easily be recovered by taking the 
limit TYij — > m u . 



1 Although we work with m d = m u , we have kept the flavor group 
to be 57/ (6 1 4) [as opposed to SU(5\3)] to allow for strong isospin 
breaking, mj ^ m u - 
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The Q field is periodic and the effects of partially 
TwBCs have been shuffled into the covariant deriva- 
tive IJ), where they have the form of a uniform U{1) 
gauge potential B^: = + iB^. The quark 

flavors are charged differently under this U(l) poten- 
tial, specifically the flavor matrix is given by = 
diag(B« ~B;, -B% Bp 0, 0, B», -B%, -B% B$. Be- 
cause the (7(1) field is uniform, the act as flavor 
dependent field momenta, each having the form B q = 
(6 q /L, 0), where L is the spatial length of the lattice and 
q is a flavor index. Notice there is no fourth compo- 
nent to B^ because the boundary conditions are only spa- 
tially twisted. The twist angles 6 can be implemented on 
the lattice by modifying the links, Uj(x) — ► Uj(x)e l6 i' L , 
where j is a spatial index. Notice that the twist angles 
vanish for the sea-sector, this corresponds to the par- 
tially twisted scenario where existing gauge configura- 
tions can be post-multiplied by the uniform U(l) gauge 
potential. We additionally must choose the ghost quarks 
to be degenerate with their valence counterparts, and 
further to satisfy the same boundary conditions as their 
valence counterparts. The determinant thus arises only 
from the sea quarks. The graded symmetry of SU(6\A) 
hence provides a way to write down a theory correspond- 
ing to partially twisted, partially quenched QCD, see [13] 
for a more rigorous discussion. 

The low-energy effective theory of partially quenched 
QCD is partially quenched chiral perturbation theory 
(PQxPT), which describes the dynamics of pseudo- 
Goldstone mesons arising from spontaneous chiral sym- 
metry breaking. In finite volume, we restrict our analysis 
to the p-regime [18| in which the long-range fluctuations 
of the Goldstone modes cannot conspire to restore chi- 
ral symmetry. The effective theory is written in terms 



of the periodic coset field E = exp(2i$//). Here the 
meson fields $ are embedded non-linearly, and in our 
conventions the parameter / is the pion decay constant, 
/ - 130 MeV. 

To obtain the relevant generalization of PQ%PT for 
Eq. |T]), we take E to transform as E — > LEi?^ un- 
der a chiral transformation (L,R) G U(6\A)l J7(6|4)a 
and write down the most general chirally invariant La- 
grangian. Furthermore, we include the uniform gauge 
potential by requiring the theory be invariant under lo- 
cal U(l)v phase rotations, see [6(. In terms of E, the 
Lagrangian of PQxPT is 

C = L- str (l^Si^Et) - Astr (m^E + EVq) + ^o$o, 

(2) 

where we have written down only the leading order terms 
in an expansion in itiq and p 2 , where p is a momentum. 
Here 3>o is the flavor singlet field $o = str(E)/\/2, and 
the action of the covariant derivative is specified by 

D^E = M E + E] + M£[T a , E]. (3) 
We have included an isovector source field Af, which will 
be utilized below. The meson modes are contained in $, 
which is a ten-by-ten matrix of fields. It has the form 

[M vv M vs xL\ 
$ = \M SV M ss x]s ■ (4) 

\ Xgv Xgs MggJ 

The mesons of M vv (M gg ) are bosonic and are formed 
from a valence (ghost) quark-antiquark pair. These ma- 
trices have the form 



V21 V22 n 21 7r 22 
"a *12 Vl2 
\^21 n 22 7 ltl V22J 



and M„ 



fvu V12 ^11 ^i 2 \ 

V21 V22 n 21 7r 22 
Kl *12 Vll Vl2 



The rjfj (vlj) mesons have quark content 77^ 



\^21 7r 22 V21 

qi<ij), while the 7r^ i^tj) mesons have quark content 



lidj (tt^j ~ Uidj ). The valence-sea (sea-sea) mesons are bosonic and contained in M vs (M ss ) as 



M, 



4>jm 4>ju 2 <f>jdi 4>]d 2 
4>lui 0/« 2 4>ldi 4>ld 2 



and M. 



r)j njl 

Klj Tjl 



Mesons contained in Xgv {Xgs) ar e built from ghost quark, valence antiquark (sea antiquark) pairs and are thus 
fermionic. These states appear as 



Xgv 





<^{ll1i2 


4>uid 1 


4>i 1 d 2 \ 
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4>uii\ 




(VU2.U2 




0u 2 d 2 


, and Xgs = 






^diui 


( t > dlU2 
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X^dzui 


( Pd2U 2 


<^d 2dl 






\^d 2j 


4d 2 J 
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In writing the PQxPT Lagrangian we have kept the fla- 
vor singlet field <E>o as a device. Expanding C to tree level, 
one finds that mesons with quark content ~ QQ have 
masses: Wqq/ = -^i ( m Q + m Q')- The flavor singlet ad- 
ditionally acquires a mass /i 2 ,. Taking the limit /^o — * oo 
(which is warranted by the strong axial anomaly), we in- 
tegrate out the flavor singlet component. The resulting 
Goldstone manifold is SU(6\A), but only the propagators 
of the flavored mesons have simple forms. The neutral 
meson propagators have both flavor connected and dis- 
connected (hairpin) terms. The flavor neutral propaga- 
tor was derived in general for the SU(M + N\N) group 



171 ] and is identical for our case because the flavor neu- 



tral states are unaffected by twisting. We will not need 
the explicit form of this propagator in order to display 
our final results. 

Further calculation at tree level shows that mesons 
with quark content QQ have kinematic momentum Pqq' 
given by Pqq> = p + Bq — Bq< . At finite volume, 
mass splittings are generated, and kinematic momenta 
are renormalized by infrared effects These effects 

were estimated to be on the order of a few percent on 
current lattices. As we shall calculate finite volume cor- 
rections to the pion form factor, which themselves are 
at the percent level, we can safely ignore mass splittings 
and momentum renormalization below. 
Pion Form Factor in Finite Volume. To calculate 
the pion form factor on the lattice using TwBCs in the 
Breit frame, two separate current insertions are used. In 
the first (second) insertion, the current strikes the up 
(anti-down) quark. Matching to this lattice setup, in our 
theory we have 2 



(4i(o)iJ>i(o)) 



e d =o 



<*&(o)|J>£(o)) 



Tr+C-jOI^K+d,)), 



= 

(5) 



with p = 9/L fixed. The currents J*' 2 correspond to the 
photon hitting the u and d quarks, respectively: 



Jl = qddu^d 2 . 



(6) 
(7) 



Here q u , qd are the light quark electric charges. The cur- 
rent insertion effectively injects momentum by changing 
flavors, from u\ to u 2 in the case of J'h, for example. 
These quark-level operators must of course be matched 
onto the effective theory. 



2 There are three other combinations of matrix elements that lead 
to the same infinite volume physics. This is due to the redun- 
dancy of the spectator quarks, which are evaluated at zero twist 
angle. We have verified that each of the four matrix elements cal- 
culated in the effective theory yield the same answer at infinite 
and finite volume. 



At leading order the matching has already been per- 
formed in Eq. ^ as we have minimally coupled the 
source field A^. As explained in the generators T a 
can be chosen purely in the valence sector. The simplest 
choice is the obvious one: (T )ij = Si 2 Sji corresponding 
to the current j\, and {T 2 )ij — 8^5 ji corresponding to 
J 2 . At next-to- leading order there is only one additional 
term needed in the effective theory. 3 This local correction 
to the current has the form 



(5J° = zagA, 



str 



W[£„st,A,S]-[A 



(8) 

This is just a rewriting of the current derived from the 
relevant term of the Gasser-Leutwyler Lagrangian [l9j |. 
As S J° will only contribute to tree level, it will not appear 
in the expression for the finite volume correction. It is 
instructive to note that this current depends only on the 
flavor non-singlet part of the generators. Any local terms 
in the effective theory arising from flavor singlet terms 
will be inaccessible using TwBCs, as these must appear 
oc str(T a ) = 0. 

Having written down the matching of the current in 
the effective theory up to next-to-leading order, we can 
now deduce the pion form factor at finite volume. As 
mentioned above, we work in the p-regime [lei ] with time 
treated as infinite in extent. Furthermore, we choose 
to calculate the form factor from the time-component 
of the current as is done in the actual lattice calcula- 
tions, e.g. 15, [3]. The loop diagrams contributing to 
the form factor are depicted in [13J. The hairpin con- 
tributions exactly cancel, which is not surprising given 
that flavor neutral mesons are electrically neutral. As a 
check on our calculation, we have verified that the infinite 
volume form factor agrees with the known SU(4\2) par- 



tially quenched expression [20L |21| . Furthermore when 



rrij = m u , the partially quenched expression turns into 
the well-known infinite volume SU(2) result [2^]. The fi- 
nite volume modification SA4(L) = Ai(L) — M.(L = oo) 
to the time-component of the current matrix element is 
given by 



SM(L) = 2E„(8)j2^ dx 



h/2 [ j,m 2 JU 



e 2 \ 



+ I 1/2 (1 - 2x)-, mj u + 4t(1 -x) — 



L 2 I 



(9) 



where we have set the overall charge q u —Qd = 1- The pion 
energy is E^ifi) — \J m 2 + 6 2 /L 2 , and the finite volume 
modification is encoded in the function Ii/ 2 (v,m 2 ) [l3T ]. 
There is no isospin breaking term in the above expression 



3 There are additional couplings that are exactly cancelled by 
wavefunction renormalization. Not surprisingly, these can be 
removed from the current by a field redefinition. 
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FIG. 1: Finite volume shift of the pion form factor in the 
Breit frame. Plotted vs. 8 is the relative change &l[G- w {Q 2 ) — 
1], where Q 2 = 6 2 x 0.025 GeV 2 . 

as it would have to be proportional to the difference of 
the initial and final state energies. Additionally there 
is a symmetry under 9 — > —9. Both of these desirable 
features are characteristics of the Breit frame kinematics. 

Numerically we can estimate the finite volume effect. 
For ease we consider the unitary point, m 2 u = m%, and 
choose to be non- vanishing in a single lattice direction. 
We consider the relative difference 

where Q 2 = A9 2 / L 2 is the momentum transfer squared, 
and the form factor G 7T (Q 2 ) has the usual definition: 
(7r(-0)|J 4 |7r(0)) = 2E 7r (9)G 7r (Q 2 ). In Fig. Q] we hold 
the lattice size L fixed at 2.5 f m, and plot the relative dif- 
ference SL[G n (Q 2 ) — I] as a function of 9 for a few values 
of the pion mass. On current lattices, the volume correc- 
tions are at the percent level and become non-negligible 
only for light pions with small twists. 
Conclusion. Above we investigate an extension of 
PQxPT which is relevant for lattice calculations using 
TwBCs in the Breit frame. For hadrons consisting of 
two light quark flavors, the appropriate theory has an 
enlarged SU(6\4) flavor group. The additional valence 
quarks are fictitious flavors differing only in their bound- 
ary conditions. This theory is described in [14j . and we 
utilize it here to determine finite volume corrections for 
lattice calculations of the pion electromagnetic form fac- 
tor. The result we derive is quite compact, maintains 
isospin symmetry, as well as a discrete rotational sym- 
metry. This is in contrast to the result derived employ- 
ing rest frame kinematics, which possessed both isospin 



breaking and cubic symmetry breaking terms from finite 
volume effects [lijj . While numerically the finite volume 
effect is demonstrated to be small for either kinematics, 
it would be interesting to investigate nucleon isovector 
quantities in the Breit frame. Results for the nucleon 
isovector magnetic moment suggested a larger than ex- 
pected volume correction arising from cubic symmetry 
breaking terms [l2^ . The Breit frame kinematics could 
mitigate such effects. 

We thank S. Simula for correspondence. This work 
is supported in part by the U.S. Dept. of Energy, Grant 
No. DE-FG02-93ER-40762 (B.C.T.), and by the Schweiz- 
erischer Nationalfonds (F.-J.J.). 



* fjjiang@itp.unibe.ch 
t rbctiburz@umd.edul 
[1] T. DeGrand and C. DeTar, Lattice Methods for Quantum 

Chromodynamics (World Scientific, 2006). 
[2] J. Zinn- Justin, Int. Ser. Monogr. Phys. 92, 1 (1996). 
[3] M. Guagnelli et al., Nucl. Phys. B664, 276 (2003). 
[4] P. F. Bedaque, Phys. Lett. B593, 82 (2004). 
[5] G. M. de Divitiis, R. Petronzio, and N. Tantalo, Phys. 

Lett. B595, 408 (2004), hep-lat/0405002. 
[6] C. T. Sachrajda and G. Villadoro, Phys. Lett. B609, 73 

(2005), hep-lat/0411033. 
[7] P. F. Bedaque and J.-W. Chen, Phys. Lett. B616, 208 

(2005), hep-lat/0412023. 
[8] B. C. Tiburzi, Phys. Lett. B617, 40 (2005). 
[9] J. M. Flynn, A. Juttner, and C. T. Sachrajda (UKQCD), 
Phys. Lett. B632, 313 (2006), hep-lat/0506016. 
[10] D. Guadagnoli, F. Mescia, and S. Simula, Phys. Rev. 

D73, 114504 (2006), hep-lat/0512020. 
[11] G. Aarts, C. Allton, J. Foley, S. Hands, and S. Kim, Nucl. 

Phys. A785, 202 (2007), hep-lat/0607012. 
[12] B. C. Tiburzi, Phys. Lett. B641, 342 (2006). 
[13] F. J. Jiang and B. C. Tiburzi, Phys. Lett. B645, 314 

(2007), hep-lat/0610103. 
[14] P. A. Boyle, J. M. Flynn, A. Juttner, C. T. Sachrajda, 

and J. M. Zanotti, JHEP 05, 016 (2007). 
[15] S. Simula (ETMC), PoS LATTICE2007, 371 (2007). 
[16] P. A. Boyle et al. (2008), 0804.3971. 
[17] S. R. Sharpe and N. Shoresh, Phys. Rev. D64, 114510 

(2001), hep-lat/0108003. 
[18] J. Gasser and H. Leutwyler, Nucl. Phys. B307, 763 
(1988). 

[19] J. Gasser and H. Leutwyler, Ann. Phys. 158, 142 (1984). 
[20] D. Arndt and B. C. Tiburzi, Phys. Rev. D68, 094501 

(2003), hep-lat/0307003. 
[21] T. B. Bunton, F. J. Jiang, and B. C. Tiburzi, Phys. Rev. 

D74, 034514 (2006), hep-lat/0607001. 
[22] J. Gasser and H. Leutwyler, Nucl. Phys. B250, 517 

(1985). 



